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We investigate the continuum limit scaling of the scalar condensate in the Nf ~ 2 Schwinger 
model on the lattice. We employ maximally twisted mass Wilson fermions and overlap fermions. 
We compute the scalar condensate by taking the trace of the propagator (direct method) and by 
utilizing the integrated Ward-Takahashi identity. While the scalar condensate comes out consis- 
tent using these two methods for a given kind of lattice fermions, we find -quite surprisingly- 
large discrepancies for the scalar condensate between twisted mass and overlap fermions. These 
discrepancies are only resolved when using the point spht current for twisted mass fermions. 
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1. Introduction 

The Schwinger model [jl|] is a good test ground for 4-dimensional QCD due to the properties of 
asymptotic freedom and the existence of non-perturbatively generated bound states. In refs. ||^ we 
have studied on a lattice the scaling properties of meson masses for a number of different fermion 
discretizations using Wilson, hypercube, Wilson twisted mass (TM) and overlap (OV) fermions. 
In these investigations, we found that the pseudo scalar mass scales with an 0{a^) behaviour for 
all fermion discretizations. The same 0{a^) scaling behaviour was also observed in [^] in a finite 
volume scaling analysis in the Schwinger model. These findings can be attributed to the super- 
renormalizability of the model. 

The aim of the present work is to perform another scaling test of the scalar condensate as an 
additional, non-trivial quantity. In the case oi Nf = 1, ref. [Q] shows that the continuum and the 
chiral limit come out to be consistent with the predictions in the continuum theories [§, ||]. Here we 
investigate the scaling toward the continuum limit in the Nf = 2 Schwinger model using maximally 
twisted mass fermions and overlap fermions. Emphasis will be put on the comparison of different 
methods and definitions of the scalar condensate used for its computation. The methods we apply 
and which will be detailed below are the direct calculation through the trace of the (inverse) Dirac 
operator employed and the integrated axial Ward-Takahashi identity. As a quite surprising outcome 
of this investigation we find that a naive approach of computing the scalar condensate by using the 
local definition, {\j/{x)y{x)), does not lead to a consistent continuum limit. 

Another aspect of our work which was presented in the poster has been the investigation of 
the question of how many eigenmodes of the used lattice Dirac operator are needed to approximate 
the pseudo scalar correlator to a certain precision. In particular, we have studied how this number 
of eigenmodes scales towards the continuum limit in the case of overlap fermions. We can, for 
lack of space, not discuss this issue in this proceedings write-up and refer to ref. [^ for a detailed 
discussion. 

2. Lattice actions and calculation methods 

We have employed maximally twisted mass Wilson fermions and overlap fermions as chirally 
improved and chirally invariant formulations of lattice fermions, respectively, in this work. Since 
both of these kind of lattice fermions are 0{a) improved, we expect only 0{a^) lattice artefacts. 

The Neuberger operator [|l] as a realization of overlap fermions is given as 



mr,a \ 1 

Doy={l-^)Do + my , Do = - 



. ^kernel 



^kernel^kemel 



(2.1) 



In the following, we will use as kernel the hypercube operator [^, i.e. Dkemei = ^hyp with param- 
eters obtained from optimizing scaling. This lattice fermion has an exact (lattice) chiral symmetry 



due to the Ginsparg-Wilson relation [10]. 



For Wilson twisted mass fermions, the lattice Dirac operator is given by 
1 2 

(2.2) 
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We use the standard Pauli-matrices a^ (/i = 1,2,3) with Os = diag(l, — 1) and denote with /t the 
unit vector shift in direction jj.. The parameters mo and ;Utm denote the untwisted and twisted bare 
fermion masses, respectively. When mo is tuned to a critical value, mo = mo,crit, by tuning the 



PCAC quark mass to zero, we reach maximally twisted mass fermions [|11|]. In this case, the theory 
is automatic 0(a)-improved and physical (parity even) quantities scale with an 0{a^) behaviour 
toward the continuum limit. 

There are two methods on the lattice, the direct and the integrated Ward-Takahashi identity 
methods, to calculate the scalar condensate which we will use for overlap and maximally twisted 
mass fermions. 

The direct method is obtained from the trace of the fermion propagator with gauge back- 
grounds. The direct method is defined as follows for the overlap and twisted mass fermions respec- 
tively; 



1 



aDo, 



1 



SdLct = ^lTr (l--^)D-; , rrct = 77lTr[/cT3T3A-n!](^_^). (2.3) 

{x,x) X 



For Wilson fermion, there is a term proportional to ^ in 2-dimensions due to the explicit 



breaking of chiral symmetry [12]. However overlap fermions and maximally twisted mass fermions 



do not have such a problem and thus no ^ term in 2-dimensions [g] appears. The same is true for 
the X term in 4-dimensions. 

a' 

Next we introduce the integrated Ward-Takahashi identity method [ [l2| | which reads 

SrwT = 2m,^(P+(x)/'-(0)) , Z:i^T = 2MtmI(P+WP-(0)>. (2.4) 

X X 

2^°WT ^^ obtained from the PCAC relation with the operator P^{x) = ^{x)a-iX±[{\ — ^)i^](;c) and 
E|^T from PCVC relation for TM fermions; d:^ {V+ {x)p- {Q)) = 2Ht^{P+ {x)p- (0)) - 5^,o(5o(;c)) 
where P'^ix) = i/ir(x)a3T±i//'(x), So{x) = i\j/{x)a3Z3Y{^)- 

3. Calculating the scalar condensate 

We have carried out numerical simulations in the following setting. The lattice size is 20 x 20 
and the statistics is over 1000 thermalized and statistically independent configurations. The quark 

mass is fixed as z = f niq^J^ I = 0.4. The error estimate is done by the method of ref. [|13|]. The 



gauge action (5g) is the Wilson plaquette action with the dimensionless coupling constant j8 = -4-^. 

To obtain results for dynamical fermions, for OV fermions, we simulated only the gauge action 

and used the re-weighting method and the spectral representation from eigenvalues and eigenmodes 



..Nf 



to compute physical observables: (^)unquench = —^ — n, \, ^° ■ For TM fermions, we also generated 

(det ./ (D/))sq 

configurations for the full action using the HMC algorithm. 

We want to remark first that all results obtained by using either Sj^^gct (^direct) '^^ ^IwT' (^Iwt) 
came out to be completely consistent when considered for each kind of lattice fermion separately. 
Therefore, we discuss in the following only one of these cases. From a dimensional analysis, we 
expect a logarithmic term in j8, 0, ^ 

y^L = A + B/p + Clog{p). (3.1) 



The two flavour Schwinger model 



Kei-ichi Nagai 



Because of universality and super-renormalisability, the coefficient multiplying the logarithmic 



term is universal and can be evaluated as C 
condensate, 



2;r 



-. This allows us to define a subtracted scalar 



v/]8Esub = V^I-Clog(j8). 



(3.2) 



Fig. m shows L^y^ and ZsubA/jS in the case of OV fermions. When the logarithmic term is sub- 
tracted, we observe a perfectly linear behaviour of 'Lsub\/^ ^& a function of l/j8. 
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Figure 1: y^H and y^Hsub in the case of OV fermions as a function of 1 //3 for a fixed value of z = 0.4. 



In order to tune to maximal twist we have determined the critical Wilson mass in the pure 
Wilson theory, i.e. setting /itm = 0. by tuning the (untwisted) PCAC quark mass to zero. The 
values of the moxnt can be found in ref. [0]. 

In fig. 0, we compare the (direct) TM fermion scalar condensate with the corresponding (di- 



rect) OV fermion data. The important result is that Atm 7^ Aqv, see eq. (3.1). Thus the continuum 
extrapolated values of the scalar condensate from both discretizations differ when the fitting func- 



tion of eq. ( p.l\ ) is used. 

In order to shed some light on the discrepancy between the results for TM and OV fermions 
discussed above, we suggest to use an improved current for TM fermions, namely the 1 -point 
splitting current, given as 



\ d I _ 



r-imp/ 



(3.3) 



for the direct method. For the integrated Ward-Takahashi identity we find 



d;{V+{x)p-{0)) = 2Atta(P+(x)p-(0)) + SpcvciP^m , (3.4) 
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Figure 2: Comparison of the scalar condensate using OV and TM fermions as a function of l/j3. In the inlet 
data on larger lattices are shown, demonstrating the smallness of finite size effects. 

where 5pcvc is the chiral rotation in the twisted mass formulation and 

P^(y) = lL^ [v{y)r5r±U^{yMy + fi) + ¥{y)r5T±Ul{y-fiMy-fi)] ■ (3.5) 

" ;U=1 ^ 
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Figure 3: The scalar condensate for TM fermions as a function of l/j3, using both, un-improved and 
improved currents. 



Fig. H shows the comparison using the improved currents in eqs. (3.3, 3.5) with the un- 



improved currents in eqs. (2.3, 2.4). We find that employing the improved current, the values 
of the scalar condensate from TM and OV fermions approach each other. The effect of the im- 
proved current is illustrated in fig. ^, which shows the correlator of the pseudo-scalar meson using 
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the un-improved and the improved currents for the case of TM fermions. It can be clearly observed 
that the short distance behaviour is significantly altered leading to the different behaviour of the 
scalar condensate. 
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Figure 4: The correlator {P{t)P{Q)) and {P{t)P{0)) for twisted mass fermions. 

In fig. ^, yjSZsub of eq. ( |3.2[ ) is plotted as a function of 1//3. From this figure, we see that 
Y^Zsub shows the expected 0{a^) scaling. In addition, the results for TM fermions using the 
improved current is very close to the one obtained from OV fermions. However, the continuum 
Mmit values of the scalar condensate are still not fully consistent. 
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Figure 5: The subtracted condensate for fixed z = (m^ y^)^ — 0-4 
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4. Summary 

We have performed numerical simulations in the 2-flavour Schwinger model thought as a test 
ground for lattice QCD. In our work, the re-weighting, the HMC algorithm for the calculation of 
the fermion determinant and the spectral representation for the correlation function are used. 

We carried out a scaling test of the scalar condensate for OV and TM fermions. Quite unex- 
pectedly, we found a large discrepancy for the values of the scalar condensate between TM and 
OV fermions, even in the continuum limit. The discrepancy is particularly significant when an 
un-improved (local) current is used for TM fermions. We have suggested to use a 1 -point split- 
ting current which indeed improves the situation but does not seem to remove the discrepancy 
completely. 

As a curious observation, we remark that we have also computed the scalar condensate by 
using an un-improved definition for OV fermions defined by P{x) = v/(x)75 va(x), in which there is 
no factor (1 — ^). In this case, the value for the scalar condensate is almost the same value as the 
one obtained from the un-improved current in the case of TM fermions. 

The reasons for the quite surprising outcomes of our investigation are presently explored. 
Clearly, if our findings are not a specialty of the 2-dimensional Schwinger model but generalize to 
lattice QCD, this would have serious implications for the determination of the scalar condensate in 
QCD. 
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